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Abstract 

In this study, nonlinear partial differential equations 
governing vibrating beam- columns with variable cross- 
section of the new method is applied to solve. Our purpose is 
to enhance the ability of solving the mentioned nonlinear 
differential equation with a simple and innovative approach 
which was named Akbari-Ganji's Method or AGM. 
Comparisons are made between AGM and Numerical 
Method (Runge-Kutte 4th). The results reveal that this 
method is very effective and simple, and can be applied for 
other nonlinear problems. It is necessary to mention that 
there are some valuable advantages in this way of solving 
nonlinear differential equations and also most of the sets of 
differential equations can be answered in this manner which 
in the other methods they have not had acceptable solutions 
up to now. We can solve equation(s) and consequently there 
is no need to utilize similarity solutions which makes the 
solution procedure a time-consuming task. According to the 
afore-mentioned assertions which will be proved in this case 
study, the process of solving nonlinear equation(s) will be 
very easy and convenient in comparison with the other 
methods. According to the explanations given about the 
capabilities of this method (AGM), non-linear partial 
differential equations governing the vibration of the beam- 
column are investigated for dynamic systems. 
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Introduction 

Many physical phenomena are modeled by nonlinear 
differential equations in order to have more 
opportunities to control the real objects in the universe, 
for example, vibration of structural systems associated 
with nonlinear properties in civil engineering is 
considered to be in this category. To better understand. 


it is better to say that nonlinear oscillator models have 
widely been used in many areas of physics and 
engineering and play an important role in mechanical 
and structural dynamics (eg, earthquakes) for the 
comprehensive understanding and accurate prediction 
of motion. In accordance with the afore-mentioned 
explanations like in this case study, attempts have 
been made in order to solve this kind of problems 
which have attracted much attention in the process of 
weaving and also to find their suitable solutions. It is 
noteworthy that many of the classic analytical methods 
such as traditional perturbation techniques which 
were presented for solving nonlinear equations not 
only have many shortcomings, but also not valid for 
strongly nonlinear differential equations. Since, a 
minority of scientific problems have precise analytical 
solutions and most of the afore-mentioned methods do 
not have the capability to solve nonlinear differential 
equations easily and accurately so in this case study in 
order to finish this difficulty, a new approach (AGM) 
for solving complicated nonlinear differential equations 
in various fields of study especially in mechanical and 
civil engineering has been presented. 

In this paper, the solution of a nonlinear partial 
differential equation of a beam-column which deforms 
due to a shear force and a rotational inertia is 
investigated. In this analysis, the effects of small ratio 
of length to depth of the beam cause a considerable 
changing on the dynamical and vibrational response of 
the beam. Figure 2 obviously depicts the effects of 
shear force and inertial resistance on rotational 
acceleration of the beam cross section. As a result, in 
condition that the ratio of depth to span is smaller, two 
factors play an important role in dynamic response of 
beam-column. The first factor is deformation which 
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due to shear force and the second factor is inertia 
resistance against the rotational acceleration of beam 
cross section. According to Fig. 2. two components of 
displacement of beam cross section are as follows: 

-Rotational inertia (mi) will be caused by rotational 
section of beam (angle a ) from it's main condition the 
perpendicular condition. 

-Deformation that due to the shear stress (by 
assumption that beam cross section also remain as a 
plate) will make angle ft . 

As it has been stated that this theory has been utilized 
in fields of mechanical and civil engineering specially 
in earthquakes analysis. Also, such beams (ie beams 
connected) can be used in structural steel and concrete. 


Nomenclature 


L 

Length of the beam - column 

E 

Young's modulus 

N 

Axial force 

m 

Mass per unit length of beam - column 

k 

Shear modulus cross beam - column (For rectangular k 
= 5/6) 

A 

Cross - section beam - column 

G 

Young's modulus, shear 

qM 

Wide load on the beam - column 

C 

Damper coefficient 

r 

Radius of gyration of the cross - section( ) 

I 

Moment of inertia 

U=U(x,t) 

Vibrating beam displacement 

* 

Shape modal 

z(t) 

Displacement function of time 

P 

shear deformation angle 

a 

Rotationlnertia angle 

nij 

Rotational inertia per unit length of the beam - column 


The analytical Method 

In general, vibrational equations and their initial 
conditions are defined for different systems as follows: 


f(ii,u,u,F 0 sin (co 0 t) = 0 

(1) 

{w(0) = A, ii(0) = 0} 

(2) 


Choosing the Answer of the Governing Equation for 
Solving Differential Equations by AGM 

In AGM, a total answer with constant coefficients is 
required in order to solve differential equations in 
various fields of study such as vibrations, structures, 
fluids and heat transfer. In vibrational systems with 
respect to the kind of vibration, it is necessary to 
choose the mentioned answer in AGM. To clarify here, 
we divide vibrational systems into two general forms: 


1) Vibrational Systems without any External Force 

Differential equations governing on this kind of 
vibrational systems are introduced in the following 
form: 

f(ii, u, u) = 0 (3) 

Now, the answer of this kind of vibrational system 
is chosen as: 

u(t) = e~ bt {Acos(cot) + B sin(cot)} (4) 

According to trigonometric relationships, Eq.(4) is 
rewritten as follows: 

u(t) = e~ bt (acos(cot + cp)) (5) 

It is notable that in the above equation a = 

\ Ia 2 +B 2 and cp - arctan (— ) . 

A 

Sometimes for increasing the precision of the 
considered answer of Eq.(3), we are able to add 
another term in the form of cosine by inspiration of 
Fourier cosine series expansion as follows: 

u(t) = e~ bt {a cos (cot + cp) + d cos (2cot + (p)} (6) 

In the above equation, we are able to omit the term 
(e~ bt ) to facilitate the computational operations in 
AGM if the system is considered without any 
damping components. 

Generally speaking in AGM, Eq.(5) or Eq.(6) is 
assumed as the answer of the vibrational 
differential equation (3) that its constant coefficients 
which are a, b, d , co ( angular frequency ) and cp 
(initial vibrational phase) can easily be obtained by 
applying the given initial conditions in Eq.(2). And 
also the above procedure will completely be 
explained through the presented example in the 
foregoing part of the paper. 

It is noteworthy that if there is no damping 
component in the vibrational system, the constant 
coefficient of b in Eq.(5) and Eq.(6) will automati- 
cally be computed zero in AGM solution procedure. 

On the contrary, the parameter b in Eq.(5) and 
Eq.(6) for the other kind of vibrational system with 
damping component is obtained as a nonzero 
parameter in AGM. 

2) Vibrational Systems with External Force 

In this step, it is assumed that the external forces 
exerted on the vibrational systems are defined as: 

F(t) = F 0 sin (co 0 t) (7) 

As a result, the differential equation governing on 
the vibrational system is expressed like Eq.(l) as 
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follows: 

/( ii, u , u , F 0 sin(&> 0 f) = 0 (8) 

The answer of the above equation is introduced as 
the sum of the particular solution (w ) and the 

harmonic solution (u h ) as follows: 

u h (t ) = e~ bt { A cos (cot) + B sin (cot)}', 

(9) 

u p (t) = M cos ( co 0 t) + N sin (co 0 t ) 

Then 

u(t) = u p +u h (10) 

By utilizing trigonometric relationships and 
substituting the yielded equations into Eq.(10), the 
desired answer will be obtained in the form of: 


u(t) = e bt {a cos (cot + <p)) + d cos (co 0 1 + (/)} (11) 


In order to increase the precision of the achieved 
equation, we are able to add another term in the 
form of cosine by inspiration of Fourier cosine 
series expansion as follows: 


u(t) - e bt {a.cos(cot + (p) + c.cos(2cot + (p)} 
+ d. cos (co 0 t + <f>) 


( 12 ) 


Finally, the exact solution of the all vibrational 
differential equations can be obtained in accor- 
dance with the following equation: 

00 

u(t) = e~ at { Yj b k cos(kcot + (p k )} + dcos(co 0 t + <f>) (13) 
k = i 


To deeply understand the above procedure, 
reading the following lines is recommended. 


Since the constant coefficient ( b ) in vibrational 
systems without damping components is always 
obtained zero, we can add the term (b.t) instead of 

(< e ~ bt ) in Eq.(12) to decrease computational operations 
in the following form: 

u(t) = (b.t) + a.cos(cot + cp) + d.cos(2cot + (p) (14) 

Based on the above explanations, by applying initial 
conditions on a system without damping component, 
the value of parameter (b) is always zero for Eq.(12) 
and Eq.(13). Therefore the role of parameter ( b ) in the 
both of Eq.(12) and Eq.(13) each of which can be 
considered as the answer of the vibrational problems is 
individually considered as a catalyst for increasing the 
precision of the assumed answer. However according 
to Eq.(16) and Eq.(20) after applying initial conditions 
on the vibrational system in both states (with external 
force and without external force) by AGM, the value of 
parameter ( b ) is computed zero because the mentioned 
system has a free vibration without any damping 
component. 


Again, it is mentioned that in order to decrease 
computational operations for systems without 
damping component and since we know that (b) in the 

term ( e ~ bt ) is zero so ( e ~ ht ) can be omitted from Eq.(ll). 
Consequently, Eq.(ll) which has been considered as 
the answer of the systems without any damping 
component can be rewritten as follows: 

u(t) = a. cos (cot + cp) + d.cos(co Q t + (/)) ( 15) 

The constant coefficients of Eq.(ll) or Eq.(12) which 
are { a , b , co , cp, c,d, (/)} will easily be computed in 
AGM by applying the initial conditions of Eq.(2). 

Application of Initial Conditions to Compute Constant 
Coefficients and Angular Frequency by AGM 

In AGM, the application of initial conditions of Eq.(2) 
is done in the two following forms: 

1) Applying the Initial Conditions on the Answer 
of Differential Equation 

In regard to the kind of vibrational system (with 
external force and without external force) which 
was completely discussed in the previous part of 
this case study, a function is chosen as the answer 
of the differential equation from Eq.(5) or Eq.(6) for 
the systems without external forces and from Eq.(ll) 
or Eq.(12) for the defined systems with external 
forces and then the initial conditions are applied on 
the selected function as follows: 

u(t) = u(IC) (16) 

It is notable that IC is the abbreviation of 
introduced initial conditions of Eq.(2). 

2) Appling the Initial Conditions on the Main 
Differential Equation and Its Derivatives 

After choosing a function as the answer of 
differential equation according to the kind of 
vibrational system, this is the best time to substitute 
the mentioned answer into the main differential 
equation instead of its dependent variable (u ) . 

Assume the general equation of the vibration such 
as Eq.(l) with time-independent parameter (t) and 
dependent function (u) as: 

/(ii, u, u, F 0 sin(&> 0 r) = 0 (17) 

Therefore, on the basis of the kind of vibrational 
system, a function as the answer of the differential 
equation such as Eq.(5) or Eq.(6) and Eq.(ll) or 
Eq.(12) are considered as follows: 

u = g(t) (18) 
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In this step, the afore-mentioned equation is substi- 
tuted into Eq.(17) instead of (u) in the following 
form: 

fit) = fig" it ), g\t), g{t ), F 0 sin(&> 0 0) (19) 

Eventually, the application of initial conditions on 
Eq.(19) and its derivatives is expressed as: 
f(IC) = f(g'{IC ), g\lC), g(IC ), ...), 

/ (/C) = / (g'(IC), g(IC ), g(/C), ...), 

/ (/C) = / (g'\lC),g\lC),g(IC),...), 


To end up, it is better to say that in AGM after 
applying the initial conditions on Eq.(18), Eq.(19) and 
also on its derivatives from Eq.(20) according to the 
order of differential equation and utilizing the own 
two given initial conditions of Eq.(2), a set of algebraic 
equations which consists of n equations with n 
unknowns is created. Therefore, the constant coeffi- 
cients ( a , b, c , d , angular frequency " CD " and initial phase) 
are easily achieved which this procedure will 
thoroughly be explained in the form of an example in 
the foregoing part of this paper. 

It is noteworthy that in Eq.(20), we are able to use the 
derivatives of fit ) with higher orders until the 
number of yielded equations is equal to that of the 
mentioned constant coefficients of the assumed answer. 

Theorical Formulation 

Consider a beam-columns with non-uniform cross- 
section according to the following figure: 



FIG. 1 THE SCHEMATIC DIAGRAM OF THE PHYSICAL MODEL. 



q(x) 



FIG. 2 BEAM - COLUMN WITH THE EFFECTS OF SHEAR 
DEFORMATION AND ROTARY INERTIA AND EQUALFORCES 

According to the Fig. 2, the balance of forces on the 
element (beam-column) the vibration differential 
equation can be obtained as follows [5,8,9]: 


a* 2 


-mr 


-iEI 


d 2 U 

dx 2 

d 4 U 


, dU .. d 2 U 
) + c— + N — -■ 
8 1 dx 2 


r < ^ 

-[q(x)-m—-] 

dt 


El 


dx 2 dt 2 kAG dx 

E I d 2 U , L 8U 2 y n 

( ^-)Jn ( ) dx = 0 

2 L 8x 2 0 dt 


8 2 r , , 8 2 U 

[q{x)-m — r ] 


dt 2 


(21) 


The beam cross-section moment of inertia (I) and 
cross-section (A) at along the beam length, we consider 
Function as the following: 

A = A 0 (l-£x) 2 ,I = I 0 (l-sx) 4 (22) 

A 0 and I 0 is the basic beam-colum characteristics and £ 
is Constant. 


Using separate equations U = <p(x). z(t) and put in a 
non-linear differential Eq. 21, proceeds obtained as 
follows [ 10,11]: 

— - — {(2 EeIqX + JcGAq-SIcGAq- E l 0 -E£ 2 I 0 x 2 )<f>"}z 
kG\(j) 

+ — Z +— {EI 0 (l-£x) 4 <f> (IV) +8EI 0 £(£ 3 x 3 -3£ 3 X 
m m(j) ^23) 

-3£ 2 X 2 + 3ex - 1)/' + (N + 12 E£ 2 I 0 +12E/ 0 £ 4 X 2 )(f>"}z 

EA^ftiWfdx^ 3 EI 0 £ 2 x 2 </> IcAqG q(x) 

-{ \z H = 0 

2mL(f) rnkGA^c/) 


And applying the Galerkinmethod[12,13], the Eq. 23 is 
obtained as follows : 

f(t)=aZ + j3Z + rZ-AZ 3 +11 = 0 (24) 

Where dot denotes differentiation with respect to time 
and a, [3, y. A, q are as follows: 

a = — ) L . 2 , [lo L (2 e Iq£x + kG\ - r 2 kG\ 

^A)Jo </> dx L (25) 

-EI 0 -E£ 2 I 0 x 2 )f'ipdxj 



k = 


El o 

2niL\' () (frdx 


{<j> Jp (^ ) 2 dx}dx 


(26) 

(28) 
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T? = D— j 0 L { EI 0 £ 2 x 2 i- k\G q(x)}tf> dx (29) 

m k GAq J 0 cp dx 

In the above equations, $ is the modal shape function 
of the beam-column [14,15] and the center of the beam- 
column subjected to the following initial conditions at 
differential equation (24) are expressed in the forms of: 

Z(0) = A,Z(0) = 0 (30) 

Solving the Differential Equation with AGM 

On the basis of the given explanations in the previous 
section, the answer of Eq.(24) in AGM is considered as 
follows: 

z(t ) = e~ at {b cos (co.t + cp)) (31) 

In AGM, the constant coefficients of Eq.(31) which are 
co ( angular frequency ) , a , b and <p (initial vibrational 
phase) can easily be computed by applying initial 
conditions according to the physical aspects of the 
problem. 

Applying Initial or Boundary Conditions in AGM 

The constant coefficients of Eq. (31) can be gained by 
applying the initial or boundary conditions in this new 
approach. 

Due to the proposed physical model, there are no 
boundary conditions so the constant coefficients of Eq. 
(31) are obtained just with respect to the given initial 
conditions which have been presented in Eq. (24). 

It is noteworthy that in the proposed method, initial or 
boundary conditions are applied in 2 ways as follows: 

1) In general, the initial conditions are applied on 
Eq. (31) in the form of: 

6 = 0(IC) (32) 

For simplicity, IC is considered as the abbreviation 
of the initial conditions. 

As a result, applying the initial conditions on Eq.(31) 
is done as: 

z(0) = A so b cos cp = A (33) 

z'( 0 ) = 0 therefore b(a cos cp + co sin cp) = 0 (34) 

2) The application of initial or boundary conditions 
on the main differential equation which in this case 
is Eq.(24) and also on its derivatives is done in the 
following general form: 

/(z(0) -»/(z(/C)) = 0, / '(z (IC)) = 0 , (35) 

Therefore, after substituting Eq.(31) which has been 
considered as the answer of the main differential 
equation into Eq.(24), the initial conditions are applied 


on the obtained equation and also on its derivatives on 
the basis of Eq.(31) as follows: 

/(z( 0 )) : b (aa 2 -aco 2 -pa + y p) cos (cp) 

(36) 

+b co(2aa - fS) sin(<p) - A b 3 cos 3 (cp) + 77 = 0 

Then for the first derivative of the obtained equation, 
we have: 

/ (z(0)) = 0 

b (3aaco 2 - aa 3 + J3a 2 - ft co 2 + y a) cos{cp) 

2 2 (37) 

-\-bco(aco -3 aa + 2a [5 - y) sin(<p) 

+3 2b 3 ( a cos cp + co sin cp) cos 2 cp = 0 

By solving a set of algebraic equations which is 
consisted of four equations with four unknowns from 
Eqs.(33, 34) and Eqs. (36, 37), the constant coefficients 
a, b, cp and co from Eq.(31) can easily be yielded as 
follows: 

To simplify, the following new variables are introduced 
as: 

y/ = 4 yAa + 4 qa - 4 A A 3 a (38) 

Therefore, the constant coefficients of Eq.(27) are 
achieved as follows: 


a = ——,b = A I ^ 2 (39) 

2 a \y/-p 2 A 

And with respect to the above procedure, co (angular 
frequency), cp (initial vibrational phase) can be 
computed in the forms of: 

1 L-/? 2 A _ lf \y/-p 2 A^ 

C ° = ~la\~A. ,( P = ts m 

After substituting the obtained values from Eq.(39) 
and Eq.(40) into Eq.(31), the solution of the mentioned 
problem will be obtained as follows: 


z(t) = A / 


I if/ - p 1 A 2a \ y/ 


V ¥ 

By selecting the physical values below: 

A = 0.1,a = 1.8,/? = 0.1,r = 0.1, 

(42) 

rj = 0.00002 , A = 0.3 v ' 

the solution of Eq.(31) and the angular frequency 
which in this case study have been obtained by Eq.(41) 
and Eq.(42) are rewritten, respectively as follows: 

® = 0.2307(^> 

sec (43) 

z(f) = 0.1 00722 e -0 ' 027778 ' cos(0.2307f -0.1 1982) 

Then, the charts of the obtained solution and its phase 
plane are drawn in the following: 
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FIG. 3 THE OBTAINED SOLUTION BY AGM 




t 


FIG. 5 THE DERIVATIVE OF THE OBTAINED 
SOLUTION BY AGM 

Numerical Solution 

With regard to the given physical values from Eq. (40) 
and the determined domain t e { 0 , 80 } which is 
defined in terms of second (sec) , the solution of the 
differential equation is expressed numerically in the 
following table: 


TABLE 1 THE OBTAINED RESULTS FOR Z (t) AND ITS DERIVATIVE 
ACCORDING TO THE GIVEN PHYSICAL VALUES. 


t 

0.0 

16 

32 

48 

64 

80 

Z(t) 

Num 

0.1000 

-0.0582502 

0.0205697 

-0.0017143 

-0.0107928 

0.01021875 

Z(t) 

AGM 

0.1000 

-0.058702 

0.0230721 

-0.010931 

-0.008300 

0.0096100 

Percent 

error 

Num, 

AGM 

0.0000 

0.00776 

0.1216 

1.637 

0.232 

0.069 

dz(t ) 
dt 

Num 

0.0000 

0.0081833 

-0.0089258 

0.0061275 

-0.0028106 

0.0004500 

dz(t ) 
dt 

AGM 

0.0000 

0.078417 

-0.0085774 

0.0061506 

-0.0032000 

0.0005010 

Percent 

error 

Num, 

AGM 

0.0000 

0.0417 

0.0397 

0.00337 

0.138 

1.156 


Comparing the Achieved Solutions by Numerical Method 
and AGM 

Due to the obtained solution from Eq. (41) by AGM, 
which is Eq.(43), and the obtained Numerical Solution 
whose results have been presented in table 1, we have 
the following comparisons: 



FIG. 6 A COMPARISON BETWEEN THE ACHIEVED SOLUTIONS 
BY AGM AND NUMERICAL METHOD. 



■ ■ Numeric ““ AGM 


FIG. 7 COMPARING THE RESULTED PHASE PLANES BY AGM 
AND NUMERICAL METHOD 



| - AGM «■ — Xun\eric | 


FIG.8 THE RESULTED CHART OF VIBRATIONAL VELOCITY BY 
AGM AND NUMERIC. 

Difference the Obtained Solutions by AGM and Numerical 
Method 

The following charts are the absolute valuedifference 
between the results on the basis of the yielded solution 
from Eq.(43) by AGM and the results of table.l by 
Numerical Method: 
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FIG. 9 DIFFERENCE THE OBTAINED SOLUTIONS BY AGM AND 
NUMERICAL METHOD. 



FIG. 10 DIFFERENCE THE FIRST DERIVATIVE OF THE 
OBTAINEDSOLUTIONS BY AGM AND NUMERICAL METHOD 


Computation of Other Damping Vibrational Parameters 

The maximum amplitude of vibrational displacement, 
velocity and acceleration equations can be computed 
as follows: 

On the basis of the given explanations in part 2, the 
answer of the introduced differential equation is 
considered as: 

z (t) =be~ at [cos {cot + cp)} (44) 

The curve of locus for the maximum points of 
vibrational amplitudes can be achieved from Eq.(44) 
when cos (cot + (p) - ±1, therefore we will have: 

z m .J.t)=be- at (45) 

In the above equation, the constant coefficients 
a and b can be obtained on the basis of the given 

explanations. By substituting t k = — in Eq.(45) which 

co 

k^{ 1,2,3,...}, it is possible to gain the maximum 
amplitudes such as z lmax , z 2max and z 3max . According to 
the given physical values, the maximum amplitudes 
can be obtained as follows: 

z*ma x =b.e a ® (46) 

Therefore, we will have: 

k = 1 — » t x - — then 

® (47) 

-a(-) -0.02777 ( ) 

^imax ~ b.e °> = 0.100722c °- 230712 =0.069 

2 71 

k - 2 —> t 2 - — as a result 

CO (48) 

-a(— ) -0.02777 ( — — — ) 

^2 max — b. c ®=0.100722e °- 230712 = 0.0472 

The damping factor in the vibrational solution of 

b{e~ at .cos(&>t+ cp)} is e~ at since in the vibrational 


differential equation the term e ^ cot is the main factor 
for damping, (£ ) can be obtained as follows: 

e ~at ^ e ~^t ^ £ = ± ( 49 ) 

CO 

In Eq(49) a and co Are evaluated by applying 
boundary conditions. 

And the initial vibrational phase can be obtained as 
follows: 

V=tg-\Z) (50) 

In Eq.(50), the parameters a and co have been gained 
in accordance with the given physical values. 
Therefore in this case study, it is clear that 
<^ = 0.060091. It is possible to simultaneously depict 
the solution of the differential equation and the chart 
of locus for maximum displacements as: 



FIG. 11 THE CHARTS OF THE OBTAINED SOLUTION AND THE 
LOCUS OF MAXIMUM DISPLACEMENT POINTS 


Furthermore, the maximum velocity z max can be acquired 
after taking the first derivative of Eq.(44) as follows: 
z (t) = b{-a.e~ at .cos {cot + cp)-e~ at .co. sin {cot + cp)} (51) 

As regards trigonometric relations, the maximum 
velocity (z max (0) from Eq.(51) can be obtained as 
follows: 

z ma x(0 = W« 2 +® 2 e~ at (52) 

In Eq.(52), t is defined as t = t k = ■ — E7 w hich 

2 co 

k <e{ 1, 2, 3,...}. As a result, the equation of maximum 
points of vibrational velocity is expressed as follows: 

_ (2k-\)n 

V = ifcmax -b^la 2 +Q0 2 e a 2 - (53) 

It is necessary to mention that a , b and co have been 
obtained with respect to the given physical values. 
Consequently, the maximum vibrational velocities can 
be computed as follows: 

J71 

k = \ — ^ Amax =0-017157 ( — ) and 

sec (54) 

J71 

k - 2^z 2max = 0.01175 ( — ) 
sec 

Then, the chart of locus for vibrational velocity and the 
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own vibrational velocity chart are illustrated as: 



FIG. 12 THE ACHIEVED RESULTS FOR THE VIBRATIONAL 
VELOCITY AND THE LOCUS OF THE MAXIMUM 
VIBRATIONAL VELOCITY POINTS 

Moreover, the desired equation for vibrational 
acceleration will be obtained by computing the second 
derivative of the obtained solution or the first 
derivative of the velocity equation. And also by 
utilizing trigonometric equations, it is possible to 
compute the locus for maximum acceleration points in 

t k = — as follows: 

CO 

The locus of maximum acceleration points 
— > z(t) = b.(a 2 +(o 2 ).e~ at (55) 

The maximum acceleration point 

z kmax =b.(a 2 +co 2 ).e 0J (56) 

In regard to the given physical values, the maximum 

vibrational accelerations in damping states can be 
obtained as follows: 

if] 

k = 1 ->z lmax = 0.00367 (— ) and for 

S (57) 

If] 

k = 2 -> z 2max = 0.002515 (— ) 

s 

Therefore, the charts of locus for maximum 

acceleration points and the obtained vibrational 
acceleration equation are depicted as follows: 



FIG. 13 THE RESULTED CHARTS OF THE VIBRATIONAL 
ACCELERATION EQUATION AND ITS RELATED LOCUS. 

Results and Discussion 

The chart of the angular frequency (co) from Eqs.(40) in 
terms of initial amplitude of vibration in accordance 
with the introduced variables from Eq.(43) is 
illustrated as follows: 



FIG. 14 THE CHART OF ANGULAR FREQUENCY IN TERMS OF 
INITIAL VIBRATIONAL AMPLITUDE. 


It is clear that the more amount the amplitude of 
vibration in the initial condition (IC) is, the more 
decreasing of the angular frequency is. 

The chart of damping ratio ( £ ) in terms of initial ampli- 
tude of vibration (A) from Eq.(49) is depicted as follows: 



FIG. 15 THE VARIATION OF DAMPING RATIO IN TERMS OF 
THE INITIAL AMPLITUDE OF VIBRATION (A) . 


With regard to the above figure, it is revealed that by 
increasing the amount of initial amplitude of vibration 
(A), the value of damping ratio will be increased. 
Therefore in this step, it is logical to indicate that 
inasmuch as the formula of damping coefficient is 
obtained by c = 2mco% , there is a direct relationship 
between initial amplitude of vibration and damping 
coefficient. 


Comparing the Charts of Phase Planes for Different 
Values of Amplitude of Vibration by AGM 

For various amounts of vibrational amplitudes, take 
for example A=0.1, 0.18, 0.26 and finally A=0.34, the 
resulted phase planes are illustrated by AGM as: 

0.03 ;= 

z'(t) - 0 . 02 -] 

-0.05 = 



A=0.1 A=0.1B ■ A=0.2ti A=0.34~| 

FIG. 16 COMPARING THE OBTAINED PHASE PLANES BY AGM 
ON THE BASIS OF INCREASING THE VALUES OF AMPLITUDE 
OF VIBRATION (A) 
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Assuming hinged-hinged to beam-column, shape 
function is defined as follows: 

</> n ( x ) = 4 sin(^A (58) 

Resulting in the displacement space-time is calculated 
as follows: 

U„ (x, t ) = Z(t ) • (J)q sin(^A , n = 1,2,3,... (59) 

Displacement curve is plotted in the form of three- 
dimensional coordinates: 



FIG. 17 DISPLACEMENT CURVE SPACE - TIME MODE n=2 
VIBRATIONAL. 

Conclusions 

In this paper, complicated nonlinear vibrational partial 
differential equations have been introduced and 
analyzed in details by Algebraic Method (AGM) and 
also the obtained results have been compared with 
Numerical Method (Runge-Kutte 4th). Then, the 
vibrational velocity and vibrational acceleration have 
successfully been achieved. Afterwards, the related 
equations of locus for vibrational velocity and 
acceleration have been gained and depicted 
completely. Eventually, the equation of damping ratio 
in terms of initial amplitude of vibration and angular 
frequency has been obtained perfectly. The above 
process has been done in order to show the ability of 
AGM for solving a broad range of differential 
equations in different fields of study particularly in 
vibrations. Consequently, it is concluded that AGM is 
a reliable and precise approach for solving 
miscellaneous differential equations. Moreover, a 
summary of the AGM excellence and benefits is 
explained as: By solving a set of algebraic equations 
with constant coefficients, we are able to obtain the 
solution of nonlinear differential equation along with 
the related angular frequency simultaneously very 
easily as applying this procedure is possible even for 
students with intermediate mathematical knowledge. 
On the other hand, it is better to say that AGM is able 
to solve linear and nonlinear differential equations 


directly in most of the situations, which means the 
final solution can be obtained without any 
dimensionless procedure. Therefore, AGM can be 
considered as a significant progress in nonlinear 
sciences. Inasmuch as the shortage of boundary 
condition(s) in this method for solving differential 
equation(s) is completely terminated, AGM is 
operational for miscellaneous nonlinear differential 
equations especially for the vibrations of civil 
engineering and Mechanical machinery and it is 
prospected that this method will be applied by 
enthusiastic young researchers in the near future. 
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